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Abstract: We consider a scenario in which the baryon asymmetry was created in the 
early universe during a cold electroweak transition. The spinodal instability of the Higgs 
field caused by a rapid change of sign of its effective mass-squared parameter induces 
tachyonic preheating. We study the development of Chern-Simons number in this transition 
by numerical lattice simulations of the SU(2)-Higgs model with an added effective CP- 
violating term. A net asymmetry is produced, and we study its dependence on the size of 
CP violation and the ratio of Higgs to W mass. 
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1. Introduction 

The baryon number asymmetry of the Universe, expressed in terms of the observed ratio 
of baryon number density (ub) to photon number density (n-y) Q 

^ = 6.5+°1 X 10-1° (1.1) 

is thought to be the result of high energy processes in the very early Universe. These pro- 
cesses need to break charge conjugation symmetry (C), the symmetry under the combined 
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charge conjugation and parity (CP) and conservation of baryon number (B). In addition 
it is necessary for the processes to take place during a period when the Universe was 
sufficiently out of thermal equilibrium |2| . 

Quite a few theories of baryogenesis have been proposed that can explain the order 
of magnitude of the asymmetry (1.1), see e.g. the review g. Most of these are based 
on physics beyond the Standard Model (SM) and because of our limited knowledge of 
this physics and of the history of the very early universe, it is hard to falsify a particular 
proposal. We consider it therefore important to search for a viable scenario of baryogenesis 
within the known physics of the Standard Model. We should also include its renormalizable 
extension that includes right-handed neutrino fields and Dirac and Majorana mass terms 
— which we shall call the Extended Standard Model (ESM) — to incorporate the neutrino 
masses. 

In the (E)SM all criteria for baryogenesis are met. Baryon number is violated through 
an anomaly Q and the weak interactions violate C and CP. For three or more genera- 
tions, CP violation in the quark sector is possible through the Cabibbo-Kobayashi-Maskawa 
(CKM) mixing matrix |^. Experimentally, three generations are observed and the mea- 
sured CP violation in Kaon and heavy-quark systems is consistent with the CKM matrix 
§. There may be more sources of CP violation in the neutrino sector, thus far hidden in 
the neutrino mixing matrix. The required non-equilibrium may be found in the decay of 
heavy neutrinos falling out of thermal equilibrium at temperatures T > 10^ GeV |^] — the 
leptogenesis scenario — or in the electroweak transition. 



Tachyonic electroweak transition 



It is a great challenge to develop a scheme in which the baryon asymmetry follows from 
the physics of the (E)SM at an energy scale of order 100 GeV. At this energy scale the 
Hubble rate is only about 10~^ eV and the required non-equilibrium dynamics is then to 
be caused by the electroweak transition. Such schemes fall under the unifying name of 
electroweak baryogenesis [ll[|. Most schemes rely on the electroweak transition being 

a sufficiently strong first-order finite-temperature phase transition. It was found later that 
the electroweak transition cannot be first order in the Standard Model |l^, 13], given the 
experimental bound on the Higgs mass mn > 114 GeV 0]. This has led to investigations 
of whether the phase transition can be first order in theories with an extended Higgs sector. 



such as the Minimal Supersymmetric Standard Model [14, 15]. 

Another possibility for the required out-of-equilibrium conditions is preheating at the 



electroweak scale after low-scale inflation, either through resonant preheating ]16, 17], or 
tachyonic preheating and the creation of topological defects ]17, In these scenarios the 
electroweak transition is not a regular finite-temperature transition but a cold transition, 
caused by a hybrid-inflation type coupling of the Higgs field to an infiaton field. Such 



electroweak transitions have been investigated in more detail ]19, 2C], and the tachyonic 



case appears most promising ]21, 23]. 



In such a tachyonic electroweak scenario for baryogenesis, the baryon number generated 
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during the transition is given by the anomaly equation 



B{t) = 3{Ncs{t)-Ncsm, (1.2) 

where A'^cs is the Chern-Simons number in the SU(2) gauge fields and B is assumed to be 
negligible (e.g. due to inflation) before the transition at time t = 0. A sufficient CP bias is 
needed to produce the observed asymmetry (|1.1D. 

CP violation 

A further hurdle to be overcome is the strength of CP violation. In the Standard Model 
this has been estimated as being of order |24 



J{yu - vcfivc - ytfivt - vufivd - Vsfivs - ybfivb - Vd? ~ lO-'^ (1.3) 

where yu, ■ ■ ■ , Ut are the fermion-Higgs Yukawa couplings and Q 

J = \lm{VijVMViVk*,)\ = (3 ± 0.3) x lO'^ (1.4) 



is the simplest rephasing-invariant combination of the CKM matrix V. In |25, IC] finite- 
temperature estimates were given with Uf ^ mj/T'^ and T of order of the electroweak 
transition, T ~ 100 GeV, which led to a magnitude ~ lO"^''. Clearly, if these estimates 
are to be taken seriously. Standard Model CP-violation is much too weak for baryogenesis, 
which is the current lore. 

CP violation is a subtle effect based on quantum mechanical interference, which can 
be destroyed at high temperatures. The above finite-temperature estimate can perhaps 
be justified by a high-temperature expansion and dimensional reasoning. However, the 
zero temperature estimate ( p,.3| ) seems less reliable to us, since it contradicts the measured 
magnitude of CP-violating observables. Factors such as J are hard to avoid in analytical 
calculations, but the product of Yukawa couplings in (|1.3|) could be compensated by pertur- 



bative energy-denominators or non-perturbative effects (see e.g. [27, 25, 29| and [pO|). We 



found indications of such compensation in a computation of the effective action obtained 



by integrating out the fermions |31|. 



One motivation for our work here is the possibility that at zero temperature the strength 
of CP violation in the SM is given by J, i.e. without the product of y's. This is another rea- 
son for considering scenarios based on a tachyonic electroweak transition. In the ESM with 
Majorana masses there is less rephasing-invariance and the corresponding J' = lm{V/jV-^)\ 
could be even larger. 

Sphaleron transitions 

It is important that the effective temperature after the transition is low enough for sphaleron 



processes to be strongly suppressed, otherwise B might diffuse to zero again (see |^, ^| 
for recent results on the magnitude of the sphaleron rate). Because a tachyonic transition 
occurs at a relatively low energy, which becomes redistributed over the many degrees of 
freedom by the efficiency of the preheating process, sphaleron wash-out is not expected to 
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be a problem, as already suggested in [16|. A recent study [^] of Higgs- and W-particle 
numbers shortly after a tachyonic electroweak transition indeed showed low effective tem- 
peratures. However, particle numbers in the important low momentum modes were found 
to be still rather large, corresponding to large chemical potentials [^4|, suggesting the 
possibility of a sizable effective sphaleron rate after the transition. 

This vi^ork 

In this paper we report on a study of the asymmetry generated in a tachyonic electroweak 
transition by lattice simulations of the SU(2)-Higgs model with effective CP violation. The 
transition is modelled by a rapid quench and we use a classical approximation that can be 



justified for such transitions |21, 22|; the quantum average is approximated by an average 
over a classical ensemble of initial conditions, which is evolved using classical equations of 
motion. 

The effective CP violation we use is given by a term 

K(/<t^Tr[F'^^F^,] (1.5) 

in the lagrangian ((p is the Higgs doublet, F^'^ the gauge field strength tensor and Ff^^ its 
dual). The coefficient k parametrizes the CP violation and has mass dimension —2. One 
may write it as 

- ^'^^^ (1.6) 



167r2M2 

where M is some mass scale and 5cp is dimensionless. 

There are two possible interpretations of the term ([L.5|). We can see it as a crude 
representation of the CP violation due to the fermions in the (E)SM, in which case we choose 



M = my/. In the second interpretation the term (1.5) could be the lowest-dimensional CP- 
violating operator resulting from integrating out heavy fermions (with masses beyond the 
electroweak scale) in a path integral, in a large class of theories beyond the Standard Model. 
In this case we choose M = 1 TeV. Notwithstanding the approximate nature of the model, 
which neglects the dynamics of the (E)SM fermions as well as the U(l) and SU(3) gauge 
fields, it is very interesting to obtain an estimate of the value of 5c p that is needed for 



reproducing the asymmetry (IJ.) via (L2), in either interpretation. 



With one exception |35|, simulations in 3+1 dimensions have not included CP vio- 



lation. In [35 1 there was a search for an asymmetry induced by a chemical potential for 



Chern-Simons number. Technical difficulties were reported and no asymmetry was seen. 



Simulations in 1+1 dimensions similar to what we consider here were performed in [36, 16[ 
and in more detail for the case of a tachyonic quench in |2^. In |2^ we showed that the 
final asymmetry is proportional to the applied C violation (in 1+1 dimensions, C violation 
plays the role of CP violation in 3+1 dimensions), with an interesting dependence on the 
ratio of Higgs to W masses. 

The structure of the rest of this paper is as follows: In section ^ we review scenarios for 
tachyonic preheating after inflation. In section ^ we describe the SU(2)-Higgs model with 
the added CP-violating term, and discuss the equations of motion. Section ^ is devoted to 
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the classical approximation and the details of the applied initial conditions. In section ^ we 
present our results. Finally, section ^ concludes with our estimate for the generated baryon 
asymmetry and the magnitude of 5cp- Technical details can be found in the Appendix. 
Preliminary accounts of this work have been given in 



2. Tachyonic electroweak preheating 

As mentioned in the Introduction, with regard to the magnitude of CP violation it is 
interesting to explore scenarios in which the electroweak transition took place at zero 
temperature. Let us assume that the effective mass parameter //gg in the Higgs potential 

V^ = Vo + fils4>U + H<t>U? (2.1) 

was at first positive and then turned negative ('tachyonic'), ending up at the current value 
—is^. We recall that n is simply related to the Higgs mass, iJ? = and assume for 

the purpose of discussion that fj, ~ 100 GeV. We now discuss two realizations for the sign 
change of //^g, firstly one that is instructive but that we will argue to be not viable, and 
secondly the more conventional one proposed in [|l^ ]. 

An elegant realization for n'^g could be a non-minimal coupling to the gravitational 
field, 

filf, = ^R-f,^ (2.2) 

with R the scalar curvature. As far as we know, there are no limits known yet on the 
parameter An indicative value is = 1/6, the conformal case, so it is natural to assume 
^ ~ 1/6. As an example, consider inflation in which R is initially large and positive, 
such that fs^ff is positive. After inflation R goes down and fi^g goes through zero when 
R ~ (100 GeV)^. The transition will take place somewhat later when the Hubble rate, 
which is generically of order i?^/^, has dropped sufficently for the Hubble damping term in 
the Higgs-field e.o.m. {3H(j)) to become subdominant, say if ~ 1 — 10 GeV. For H k, 1 GeV, 
the energy density is still very large, p = SH^rup (lO^GeV)^ (we use nip = (SvrG)"^/^), 
far above the electroweak scale of about (100 GeV)^. After the transition the inflaton still 
dominates the energy and it should have released this into the SM degrees of freedom before 
big bang nucleosynthesis (BBN) can take over. 

In the conventional scenario (see e.g. 1^^, section 8.3), the universe expands in matter- 
dominated fashion when the inflaton (a) has fallen out of slow roll and is oscillating in the 
minimum of its potential, while decaying into the SM degrees of freedom ('radiation'). The 
inflaton decay rate Fu should be sufficiently slow, otherwise the maximum temperature in 
radiation pO| , 

T^ax ^ g:'^'Ml^^{r„mp)'/\ (2.3) 

might rise above the electroweak transition temperature Tc ~ 100 GeV, and the sphaleron 
rate would wash out the generated asymmetry. Here (7* is the usual effective number of 
d.o.f and M^^q is the energy scale of inflation, which by the previous argument is larger 
than ~ 10^ GeV. Requiring Tmax 

< 100 GeV leads to < lO'^^ GeV. Such an extremely 
small decay rate is perhaps not unnatural if the inflaton resides in gravitational degrees 



37, 38, m 31]. 
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of freedom, and — 1/6)^ is sufficiently small. However, it also leads to a reheating 
temperature 

Trh^5*"'^'(r.mp)^/2 (2.4) 

that is smaller than 10^^ MeV, in conflict with BBN. Moreover, the entropy in a comoving 
volume produced by the decaying inflaton, S = sa^ {s is the entropy density and a the 
scale factor), increases oc a^^/^ [40|, whereas the corresponding baryon number B = ubcl^ 



is essentially conserved for T < 100 GeV. Then the crucial ratio n^/s suffers dilution by 
a huge factor {ai/a2Y^^^ = (^1/^2)^''^ ~ 10"^"^, where ti « mp/M^g is the time at the end 
of inflation and t2 ~ 1 /'^a the time the inflaton has decayed. 

The numbers look better, but not good enough, if the inflaton potential does not 
have quadratic minimum but simply falls away such that the kinetic energy dominates 



('kinetion'), as in quintessential inflation |41]. In this case the energy density falls faster 



than the radiation in SM particles, oc a . However, going through a similar analysis 
leads to practically the same Trh, an entropy growth S oc a^/^ and a dilution factor « 10~^. 
In case is negligible, the entropy would be conserved and ub/s would be constant. 
However, it then takes too long for the energy density in SM d.o.f., /jsm oc a~^, to become 
comparable to p^j. This would happen at a temperature much too low for BBN: T/Ti = 
ai/a = {psM/Pa)T - 10"'^ or T « 10"^ MeV. 

These problems are avoided if the energy scales of inflation and radiation are compa- 
rable, specifically low-scale inflation with Minfi ~ 100 GeV. For this the realization 



does not work (unless ^ > 10^^!). In |16, 17] a hybrid inflation mechanism was proposed 



in which the Higgs held is coupled to the inflaton, 

/^eflf = - (2.5) 

and a was assumed to roll from large values towards (7 = 0. Considering the effect of 
radiative corrections of the Higgs to the inflaton led jl^ to the conclusion that inverted 
hybrid inflation is a better option (see also p2|). 



/^eff = M</, - A<x0 , (2.6) 

where now a runs from to [{p? + p'^ / ^a4iY^'^ ■ An explicit example of the inflaton-Higgs 
potential is given in [p^: 

y(a, (t>) = Vo-\ Asa^ + ^Acea^ + ^-<t> ^')'^' + J^'^' (2-7) 

(for simplicity for one real (p), with A5 = 7.3 x 10"^ GeV"\ kq = 2.4 x 10"'^GeV"^, 
pI = 1000 GeV^ A^^ = 0.04, and Vq such that V = 0m its minimum, Vq ~ (86 GeV)^. 

Choosing A = 1/9 leads to the vacuum expectation values {(p) = 237 GeV (reasonably 
close to the correct 246 GeV), (a) = 426 GeV, and p = 79 GeV. The inflaton-Higgs cou- 
pling A(j(^ is quite strong, which causes substantial mixing in the minimum of the potential. 
The eigenvalues of the mass matrix are given by (74GeV)^ and (147 GeV)^, differing con- 
siderably from the diagonal elements d'^V/da'^ = (121 GeV)^ and d'^V/dcjP' = (lllGeV)^, 
so one may wonder if the SM Higgs physics is not too much affected. 



-6- 



Another question that needs further investigation is the fact that in this model the 
inflaton decay width is much larger than the Hubble rate. At this low-scale inflation 
the Hubble rate is only H ~ 10~^ eV, whereas one expects To- ~ 1 GeV. In the regime 
To- >> H, warm inflation has been advocated |43], which puts into question our assumption 



of a cold universe after inflation. 

However, the problems with the previous realization (^]^) are avoided. After the 
transition the inflaton decays rapidly into the SM d.o.f. via the mixing with the Higgs. 
Neglecting the tiny Hubble rate and using energy conservation, the 'reheat' temperature 
is approximately given by 

3U \ ' 1/4 , , 

. 2.8 



With = 86.25 for the effective number of d.o.f. below the W mass (leptons, quarks, 
gluons, photons), and a Higgs mass of, say, 160 GeV, this gives T ~ 51 GeV, well above 
the BBN and QCD transition temperatures. 

Because of the uncertainties in the specific realization of the assumed tachyonic elec- 
troweak transition, we model the transition by an instantaneous quench: 

nlff = t < 0, 

= -n^, t > 0. (2.9) 

This is a limiting case of the transition in a finite time treated in 23 1. In this way 
we 'shield' the SM from the uncertainties of the inflaton and presumably, the quench will 
generate a maximal baryon asymmetry for a given amount of CP violation. Requiring zero 
vacuum energy gives 

Vo = //4A, (2.10) 

with a 'reheat' temperature similar to the model (p. 7]). 

The initial state for the Higgs field is {(p) = with jj.'^q > 0. When /i^g < the 
Higgs field suffers the spinodal instability: (0) does not change, but its low momentum 
modes grow exponentially (2^. Soon the quartic term in the Higgs potential kicks in 
and eventually with the other SM couplings the fields will thermalize in the broken phase 



minimum. Initially there is rapid effective-thermalization 1 34 1 , called 'tachyonic preheating' 



3. The SU(2)-Higgs model 

To study the baryon asymmetry emerging after the electroweak transition we used the 
SU(2) Higgs model with effective CP violation, given by the action 



S = - (Tx 



(3.1) 



Here F^^ = d^A^ - d^A^ - i[A^,A^], D^cj) = {d^ - iA^)(j), with the Higgs doublet, 
F^u = e^upaFP''/2, our metric is (-1,1,1,1) and eoi23 = +1- As usual A^ = A^t"- 12 
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with r°, a = 1,2,3, the Pauh matrices. The vacuum expectation value of the Higgs field 
is \{4>)\ = v/\^, with V = /i/\/A, and the Higgs and W masses are given by mu = 
V, niw = gv/2; Vq = //^/4A. The magnitude of the effective CP violation is 
parametrized by n. For later reference, we define the dimensionless 

k = le-K^Km^. (3.2) 

The equation of motion for the Higgs doublet is 

D^D^' + 2A0V - KTr Pf'^F^^^ (/) = 0. (3.3) 

For the gauge fields the equations of motion are given by 

Do (^1 E^^ - 2k4>HBI^ - ekimDi (^^ + 2k4>UE^^ + jl = 0, 



(3.4) 



where E'^ = F^q and i?^ = efcimF^'^/2 are the SU(2) electric and magnetic fields, is the 
covariant derivative in the adjoint representation, e.g. 

D^Bl = d^Bl + ea,cAlBl, (3.5) 

and is the Higgs contribution to the SU(2) current 



= i{D^4>)^^—c^ - i(l>^^D^(^. (3.6) 



In addition, the Gauss constraints have to be satisfied, 

Dk E^, - 2^cl)UBt^ + Jo" = 0. (3.7) 

If these conditions hold at one time, they hold at all times as a consequence of the equations 
of motion (|3^. 

If (/)^(/) were constant, the CP-violating terms are ineffective, since then the K-term in 
the action is the integral of a total derivative, 

TtF^.F^''' = l6^d^j^s^ (3.8) 

where j'^g is the Chern-Simons current 

. e^''^^ (^AlF-, - ^e^.^AlA'lAl^ . (3.9) 



32tt 

Similarly, for constant (p'^cp, the k terms drop out of the field equations because of the 
Bianchi identities e'^^^'^ D\F^j^ = 0. Hence, an alternative version of ( |3.4D and ( |3.7| ) is 



= DoEl - 2g\B^do{(^U) - ekimDiB^ - 2g^ neuAd^ <t^) + 9^fk, (3-10) 
= DkE^^ - 2g\B^M(^U) + g^j^, (3.11) 



respectively. 
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4. Initial conditions and classical approximation 



We consider an initial state where the Universe is at zero temperature, with the Higgs field 
expectation value at zero. In the vacuum, there are quantum fluctuations, and we would 
like to use those to seed the Higgs symmetry breaking. Several ways have been suggested 
(see for instance jl^, ^). We follow the line described in |Q (see also |2l| for the more 
realistic case of a rapid but not instantaneous quench), by solving the quantum evolution 
in the limit of zero coupling, A = 0. We also neglect the interactions with the gauge fields 
for the moment. In this limit, the Higgs potential after the quench at t = is just an 
inverted parabola 

]-fi'^ct>a(t>a (4.1) 



where the (pa, ct = 1,2,3,4, are four real fields representing the complex Higgs doublet. 
For the moment we consider just one of those real components. It is straightforward to 
solve the operator equations of motion for t>Q with the initial condition that the field is 
free with mass ^ for times t < 0. In terms of Fourier modes 



-ikx 



X 



(4.2) 



in a periodic volume L^, one then finds that modes with k < ^ are unstable and grow 
exponentially. For ji^ — k'^t ^ 1, 



9k oc exp 



kH 



(4.3) 



and the particle numbers of the unstable modes in the initial state |0) (the vacuum for 
t < 0) also grow exponentially. So these dominating modes become classical and indeed, 
the expectation value of generic products of field operators can be reproduced by a classical 
gaussian distribution [|2], 




+ 



14" P 



2 i,'^ V'^fc + l/2 + nfc Hk + 1/2 - fik 

|k|<M 



(4.4) 
(4.5) 



where vr is the canonical conjugate to (p. The particle numbers Uk, fik and frequencies cok 
are defined in terms of the two-point functions 



Uk + 



1\ 1 



2 J ujk 



(TTk'/'k) = (^k'/'k) + i = nk + 



(4.6) 
(4.7) 
(4.8) 
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which can be calculated to be |22| 



Mi) 



1 + 



+2 



UJi- 



UJ 



+2 



-2 



1 sin^(a;^ t) 



1 cos^(a;^ t) 



UJ,. 



1 sm(2w. + 7T> 



(4.9) 

(4.10) 

(4.11) 
(4.12) 



Explicit expressions for n^, and are easily obtained, which show that for \uj^ \ t > 1, 
nfc + 1/2 + nfc grows exponentially, whereas the difference + 1/2 — rapidly approaches 



zero (see figure 2 in |Q). This means that the typical ^j^ grows exponentially and the 
typical — > 0. In terms of and TTk the distribution gets squeezed. 

It therefore makes sense to treat the dynamics classically, as soon as the particle 
numbers are large, ^ 1/2. This is the case for sufficiently long "roll-off" times, keeping 
in mind that the quadratic approximation will break down when the interaction terms 
becomes non-negligible. For a given choice of "roll-off" time and thus n^, and cj^, we 
reproduce generic quantum correlators of the initial conditions by an ensemble of random 
classical initial field configurations. These can then be evolved using the classical equations 
of motion. In this classical approximation the quantum averages are replaced by averages 
over the initial ensemble. 

The classical evolution can be carried out on a computer, including fully the non-linear 
interactions. This assumes of course that the particle numbers in the gauge fields will also 
grow large, such that they can also be treated classically. A large growth of gauge field 



modes could be seen a posteriori. In [34|, we found that the Higgs occupation numbers for 
the unstable modes indeed become very large (no ~ 100) during the instability. In addition, 
the gauge fields also acquire large occupation numbers, thereby supporting the classical 
approximation for the whole system, at least for the time scales under consideration here. 
The large particle numbers in the gauge field modes suggest that our major observable to 
be computed, (A'^cs(*) — -^cs(O)), can also be obtained reasonably accurately in the classical 
approximation . 

The Gauss constraints ( [j.Tl ) need to be imposed on the initial conditions. We do this 
as follows. As long as the free-field approximation is valid, the gauge fields cannot feel 
the exponential growth of the Higgs fields, so we initialize A'^ = 0. Then the covariant 
derivative in (|3.7D reduces to an ordinary derivative and it is straightforward to solve for 
the electric field components, S^, given the Higgs charge densities j^'^ drawn from the 
distribution (4.4). However, before doing so one constraint has to be imposed first: in 
finite volume the global charge vanishes, 



J d'xdkFSJg^ = 0. 



(4.13) 
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We thus modify the distribution (|4.4| ), 



PiO^PiO^iG), (4.14) 

where <5(G) encodes this global Gauss constraint. This makes the distribution not quite 
Gaussian. For details see appendix 

We will apply two types of initial conditions, dubbed 'Just-a-half and 'Thermal'. For 



additional details on these schemes, see [22 
4.1 Just-the-half distribution 

In the quadratic approximation, classical and quantum evolution is identical, and we can 



therefore choose to sample the initial distribution ( [4.4| , |4.14| ) at "roll-off" time zero, when 
the n's and w's are simply given by 



nk = nk = 0, ujk = vV + A?2. (4^5) 

leaving just the 1/2 in the denominators in ([4.4|). We only initialize momentum modes that 
are unstable and that will acquire large occupation numbers under the spinodal instability, 



so |k| < /i. This also avoids initializing modes close to the lattice cut-off scale [46|. We 
name these initial conditions: Just-the-half. 

4.2 Thermal distribution 

A different condition is obtained from an initial state that is thermal but at a low temper- 
ature T <^ niH ■ We consider this possibility to get some indication of the sensitivity of the 
results to the initial conditions. The initial particle numbers of the Higgs fields are chosen 
to be Bose-Einstein distributed: 



rifc = (exp (wfc/T) - 1) S "-fc = 0, iVk = vV + ^> (4-16) 

We use T/niH = 1/10. In this case we initialize all the modes, not only the unstable ones, 
but of course, only the unstable ones will grow. In practice the large-momentum tail is so 
suppressed that it should make little difference to introduce the cut-off k < fi. We name 
these initial conditions: Thermal. 



5. Numerical simulation 

For the numerical simulation we translated the action to a lattice in (Lorentzian) space- 
time according to the usual method of lattice gauge theory, from which the discretized field 
equations follow in the standard fashion. The Gauss constraint ( |3.7| ) is compatible with the 
discretization and since the initial Aq = 0, it stays zero and we are using temporal gauge. 
Because of the CP-violating terms the numerical algorithm for the equations of motion 
turns out to be 'implicit', causing it to be rather computer-expensive. Lattice details are 
given in the appendix. 

We evolved ensembles of 45 initial conditions in time using classical dynamics. In 
two cases (Thermal initial conditions, k = IGvr^Km^ = and 3) we averaged over larger 
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ensembles of about 300 configurations. We used volumes {Luih)^ = 21^ with periodic 
boundary conditions, and a spatial lattice spacing aniH = 0.35, so with 60'^ lattices. This is 
a compromise between computer time and capacity and the need to have enough unstable 
modes on the lattice. Here, we have about 50 such modes. Also, we found that with 
this choice of lattice spacing discretization errors in the Chern-Simons number and Higgs- 
winding number (see below) were small enough for our purpose. 



(/>t(/> = -J^ / d^xct)^, (5.1) 



5.1 Observables 

The spatial average of 

rh]rh = 

L3 

is a good indicator for the development of the instability in time. Our most important 
observable is the Chern-Simons number 

iVcs = j d'xfcs, (5.2) 

since the change in time of its expectation value determines the baryon asymmetry through 
eq. (pD, 



B{t) = 3{Ncs{t) - Ncsm = 3^* dx^ J d'x {d^J^^) = ^ dx"^ j Sx {ElBt). (5.3) 



Another interesting observable is the winding number in the Higgs field 



iVw = ^ y d^xeu^Ti {dkVV^diVV^dmVV^) , (5.4) 

V = ^^I^ E SU(2). (5.5) 

For configurations near equilibrium with low energy in the covariant Higgs derivatives, the 
gauge field is not far from being pure-gauge, A^^ ~ —id^VV'^, and then A^cs — Aw- 

Observables such as the Chern-Simons number and the winding number are notoriously 
difficult to implement on the lattice in the quantum theory, but they turn out to be man- 
ageable in our classical approximation because the high-momentum modes are suppressed 
in the initial conditions. 

5.2 Single trajectories 

We evolved the initial configurations to time t = 100 m]^^, and computed the volume- 
averaged Higgs field (j)'^4>, the Chern-Simons number A"cs and the Higgs winding number 
Aw- The initial Chern-Simons number was set to zero, so we really computed its change 



through eq. (5^), which is gauge invariant. 

Fig. H shows the observables versus time from a single such trajectory in configuration 
space. We see that "falls off the hill into its broken phase minimum" , where it performs 
a damped oscillation as the energy becomes slowly distributed over the higher momentum 
modes. Meanwhile, the Chern-Simons number and winding number bounce around, until 
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50 



100 



Figure 1: Volume-averaged Higgs field 20l0/u^, Chern- Simons number and winding number, for 
a typical trajectory with Thermal initial conditions, mn/mw — 1, and no CP violation, fc = 0. 




40 60 

Figure 2: Example of iVcs and iV^ for two trajectories with the same initial conditions, without 
(fc = 0) and with (fc — 3) CP violation. Thermal initial conditions, mn/mw — 1- 



they settle near the same integer value (for A'„ the difference with the integer is a finite- 
lattice spacing effect). The winding number appears to settle first and then the Chern- 
Simons number approaches it somewhat later. 

In Fig. ^ we show a particular initial configuration evolved with and without CP 
violation. In this case, the final Chern-Simons number and winding number are shifted 
approximately by an integer (1) under the influence of CP violation. It is these shifts that 
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will give us a net Chern-Simons number asymmetry. Due to the chaotic nature of the 
equations of motion the shifts cannot be predicted from the magnitude of the CP violation 
(for example, in figure |^ the final A'^cs — 1 for k = and ~ for k = 3). However, on the 
average, more shifts happen to one side than to the other. 




Ncs, Nw 



Figure 3: The distribution of final A'^cs and Ny^; Thermal initial conditions, niH /mw = 1, no CP 
violation (fc = 0), 303 configurations. 




Ncs, Nw 

Figure 4: As in figure ^ with CP violation, k ~ i. 

Fig. ^ shows the distribution of final Chern-Simons number and winding number. 
Although both observables cluster around integers, the winding number is more peaked. 
The relatively short evolution time {mnt^naX = 100) does not allow Chern-Simons numbers 
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to settle completely. Furthermore, the system has a non-zero effective temperature, so the 
Chern-Simons numbers need not be integer. We have checked for a few configurations, that 
the Chern-Simons number settles eventually, sometimes as late as mnt = 500. For most 
cases the Chern-Simons number is stuck already at time rriHt = 100. 

Fig. ^ shows the final distributions when adding CP violation. Note that the initial 
configurations in Figs. |3| and ^ are exactly the same. It is not obvious that anything has 
changed and to spot the asymmetry we have to form ensemble-averaged quantities. 

5.3 Ensemble averages 
1.5 

1 



0.5 





-0.5 

20 40 60 80 100 

Figure 5: Averaged observables: 2{(j)'^(j))/v^ , {Ncs) and (A^w) versus time for rriH/mw = 1, Just- 
a-half initial conditions and k — 8. 

In Fig. ||we plot the ensemble averaged Higgs expectation value, Chern-Simons number 
and winding number vs. time. We see that the average Higgs field damps more strongly. 
The Chern-Simons number oscillates with a large amplitude before settling down. The 
oscillations seem to be driven by the oscillations in the Higgs field. At the end, both 
observables settle at non-zero values. Fig. ^ shows the average Chern-Simons numbers for 
different values of k. Notice that the oscillations have larger amplitude with larger k. 

Fig. 1^ sums up our results and the dependencies on the parameters. We keep mw 
fixed and plot {Ncs)/{Lm]y)^ vs. k at mnt = 100. The dependence on initial conditions 
appears to be weak. There is a hint of a dependence on mass ratio mn/'mw- 

The dependence on k is not linear for the Thermal, niH/Tnyy = 1, case. We have 
observed similar non-linear behavior for large k (large C-violation) in the 1-|-1 D Abelian 
Higgs model, but were able to establish linearity at smaller k. Here a linear regime presum- 
ably also exists at small k, but it is unclear where it ends. We lack sufficient statistics to 
draw a conclusion. Still, a linear fit through the origin is consistent with the data for A; < 10 
(dashed lines), for the niH = ^/2mw case (x^/d.o.i = 0.25) and less so for mn = TUyy 
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Figure 6: Expectation values 2{(j>^ (j>) /v'^ and (-/Vcs) for different k; niH/mw = 1 and Thermal 
initial conditions. 
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Figure 7: Results for the Chern-Simons density all fc, nfiH/mw and initial conditions. The full 
lines represent linear fits through the two points with lowest k: fc = 0, 3 {mn = mw) and fc = 0, f .5 
{niH = \/2mw)- The dashed lines are linear fits through the origin, ignoring the data at A; = 0, 
and also at fc = f 6 for the case mn — mw- 



(xVd.o.f = 2.4), 



(iVcs) 



(0.32 ± 0.15) X 10"^ /c, mH = mw, 



(5.6) 
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= (0.79 ± 0.36) X 10"^ A;, mn = V2mw 



(5.7) 



This fit takes into account that with infinite statistics the result at k = should be zero. 
However, the large iii t^is tuh = my/ case suggests contamination from non-linear 
behavior. 

To diminish the uncertainty of where the linear regime ends we focus on the lowest 
fc-values, A: = 0, 3, niH = mw, and fc = 0, 1.5, mn = V^mw- We expect these points to 
be close to, or in, the linear regime that is of physical relevance. Furthermore, since we 
used exactly the same initial configurations for the data at /c = and at k = 1.5,3, the 
effect of CP violation is the difference between the zero and non-zero k results. An overall 
shift due to finite statistics cancels out in this difference. These fits are represented by the 
continuous lines in figure 0. The fitted slopes are given by 



: (0.79 ±0.31) X 10"^ /c, mH = mw, (5.8) 

~ 1.4x10"^ A;, mH = V2mw (5.9) 

The error in the mn = mw case is obtained by shifting the data upward with the negative of 
the k = value and combining the errors at A; = 3 in quadrature. For the case rriH = V^mu 
this would lead to a too large, not meaningful error, which we have left out. We consider as 
our best estimate for the slope the result (5.7) from the previous fit (upper dashed line in 



figure 0), keeping in mind that this is perhaps an under-estimate because of the apparent 
curvature in the data. 



5.4 Initial rise 

It is clearly difficult to predict the final averaged Chern-Simons number analytically. For 
early times, specifically during the very first roll-off of the Higgs field, we can however 
understand the asymmetry generated by the driving force of the CP-violating term, through 
an analysis similar to that used in the analog 1+1 D Abelian Higgs model |2^]. Consider 
a spatial region in which the fields are gauge-equivalent to being homogeneous and choose 
a gauge in which they actually are homogeneous. For homogeneous fields, the equations of 
motion for the gauge field read: 

d^At = A\A\A^ - AIA\A\ - \ g^A^H - 2g\Bt dt{<i>U) (5.10) 

In the Abelian Higgs model in 1+1 dimensions [p^ j, the simplicity of the C-bias term 
ensures that the gauge field does not enter into the driving force (the K-term). In the 
present case, however, its effect is suppressed by 

Bl = ]^eumeabcA\A'i^. (5.11) 

To estimate the asymmetry, we assume that the Higgs field, which may be viewed as an 
0(4) vector, rolls off the quadratic potential, say in the direction 03 : 



dU^ = t?^^^(t^H = ^le^^\ 



(5.12) 



Assume that we have solved the system for n = 0. From the simulations, we know that 
the gauge field also grows exponentially, say 

At,o = Cle'^'- (5.13) 
We perturb with the k term, and write 

AI = AI^^^kAI^. (5.14) 
Then the equation for the perturbation ^ is 

d}Al^^ = -Ht\Al, - 2g'Bl,dti^^^), (5.15) 

with 

= i^'^ki - AlA'i)6at - AlAl6ki + 2AIA\ - AfAl + i g'c^USkiSab, (5.16) 



evaluated at zeroth order in k. For small times H is dominated by the (j)'cl) term, which acts 
like an effective-frequency term for the perturbation, that will turn the initial rise induced 
by the source term oc dt{4)^(t)) into the small initial bump in figures ^ and|6|. Assuming this 
effective mass term can be neglected at early times we can integrate ( ^.151) using (|5.12| ) 
and ( |5.13|) and express the result as 

For homogeneous fields the Chern-Simons density is given by (cf. (|3.9| )) 

ft I 
ncs(t) = dx^ dofcs = ^kMc AlA\Al, (5.18) 

evaluated at time t. Averaging over initial conditions the zeroth order contribution in k 
will vanish and we get 

(ncs) = {Bl,Al,) = (5.19) 

We now assume this expression, which was derived for our homogenous patch, to be inde- 
pendent of where the patch was located, and replace {B'^(fi'^(j)) — > {B'^) {(^'^ (j)) , where (-B^) 
and {(f)^(p) stand for the average over the volume as well as over initial conditions. This 
gives for the density in W-mass units, 



(ncs) V2k {B^){ct>^ct)){2/v^ 



(5.20) 



'^w (87r2)2(l +c)2 m^niH 

Using the computed values of {B"^) and and inferring c from the measured behavior 

(i?2) tx exp(4c/it), we can compare ( 5.20 ) to the data. For c we find c ~ 0.67. 

Fig. P shows the initial bump for niH/my/ = \/2, for varying k. The inset shows a 
comparison with the estimate ( ^.20D , evaluated at a time for which {4)'^(j)) = f^/6 , the value 
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Figure 8: The initial rise for various k's. The inset shows the size of the bump versus k for the 



estimate (5.20) (hne), the results of the simulations (dots) and results shifted by the negative of 
the value at zero k (see main text); mH/rnw = \/2. 



at which the curvature in the Higgs potential changes sign and the quadratic approximation 
(5.12) will be become more strongly affected by non-linearities. Because of finite statistics, 



there is an overall shift of the results that can be read off from the A; = simulation, which 
should give zero with an infinite number of initial configurations. In this case the A; = and 
1.5 simulations started from the exact same initial conditions, and so we shift the A; = 
result to zero, and the k = 2> result by the same amount (as we did in section ^]^). These 
are the crosses in the inset. The agreement of ( 5.20| ) with the data is surprisingly good, 
given the approximations made. 

5.5 Final temperature and sphaleron wash-out 

In order for baryogenesis to be successful, we need the final temperature after symmetry 
breaking to be low enough so that an asymmetry does not get washed out by equilibrium 
sphaleron processes. The usual requirement for avoiding subsequent wash-out, namely, that 
the sphaleron energy is much larger than the temperature, can be conservatively formulated 



as ||l^ v/T > 1.49, or T/v = 0.67. The estimate (|2.8D , based on distributing the Higgs 
potential energy (Vq = over g^: relativistic degrees of freedom, leads to 

- = 0.442 — , 5.21 

V \ V / 

where we used 5* = 10 for the Higss and W degrees of freedom. With v = 246 GeV and 
a Higgs mass in the range 114 - 200 GeV the requirement is amply satisfied. For 
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example, for mn = V^fnw — 114 GeV the estimated temperature is T = 74GeV. 

In fact, our system is not yet in equilibrium after the transition and our particles are not 
massless. In |^4| we computed Higgs and W particle numbers using the field configurations 
produced in our numerical simulation, from which we obtained effective temperatures and 
chemical potentials by comparing to a Bose-Einstein distribution. We found 

rp ch 

— -0.4, -^-^ ~ 1, mH = V2mw (5.22) 

At first glance this temperature should keep us safe from equilibrium sphaleron wash-out. 
The implied temperature 0.4 x 114 = 45 GeV is lower than the 74 GeV found above, which 
can be ascribed to the occurrence of the chemical potential. However, this large chemical 
potential implies that the occupation numbers are still very large for the low momentum 
modes, = [exp(\/m^~+~A?/T — pL^^/T) — 1]~^ ^> 1. Since these modes are the ones 
responsible for sphaleron processes, the actual rate may be larger than one would expect 
from the temperature alone. 

In the Standard Model, the W and Higgs particles decay into the lighter particles, 
which will rapidly lead to a lowering of the temperature. For the W, the decay width is 
about 3 GeV and the Higgs width is expected to be similar, so that in a time span of a 
few hundred rn^ these particles have decayed. We do not see any wash-out on time scales 
of a few hundred mj^. Then the final temperature will be determined by light degrees of 



freedom and thus close to the estimate (2^) with the somewhat more favorable = 86.25: 



^ = 0.258 

V \ V 



(and a temperature T = 43 GeV for mn = 114 GeV). We do not see a problem with 
wash-out. 

5.6 Effective sphaleron rate 

In Ijl^ an effective sphaleron rate was introduced and used in estimating the effect of the 
CP-bias term using quasi-equilibrium concepts, 

where 5cp/M^ = 167r2K/3 and T^s 

is an effective temperature estimated to be about 350 



GeV. This effective sphaleron rate was computed recently |23], with parameters different 
from ours (and with a finite quenching rate), and it is interesting to see its behavior in our 
simulation. The effective rate per unit volume can be defined as 

Teff = ;^ I [{Ncsitf) - {Ncsit))'] , (5.25) 

where we set A^cs = at t = 0. In Fig. |9| we show Fefr versus time. It turns out to be 
insensitive to the choice of initial conditions, and insensitive to rriH/rnw- In equilibrium 
this quantity is the diffusion rate of Chern-Simons number, in which case it is the slope of 
a straight line approximating {N^a,) — (-Acs)^ (see inset). 
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Figure 9: The effective sphaleron rate Tcs/rnj^; mH/mw 
Inset: integrated rate dt' Vcs{t') /m^fj. 



1, /c = 0, Thermal initial conditions. 



In the tachyonic transition such a diffusion regime is preceded by a period of preheating, 
where the gauge fields acquire energy and as a result a non-zero (-/V^g). The early-time 
rate at which this happens does not have anything to do with the equilibrium sphaleron 
rate. Rather, it is determined by the time scale of the preheating mechanism. Note how 
this rate is initially in phase with the Higgs field, but it has a large magnitude for small and 
large values of (</>^</'). At these early times the system is not experiencing potential barriers 
proportional to {4>^4>), with sphalerons mediating Chern-Simons number change. Even 
at the latest time the rate is not as small as the temperature would suggest. Averaging 
out the oscillations in the region 10 < tniH ^ 40, the effective slope (rate) is about 
Tes/iTij^ 0.35 X 10^^, about the same value as obtained in [|^]. By ( ^.24 ) this leads to 
the estimate (for mn = rriy/) 



ns/m^r « 0.1 X 10 ^ k, 



(5.26) 



an order of magnitude smaller than our result (|57 
get HB/m^). 



(which should be multiplied by 3 to 



6. Conclusion 



Using numerical simulations that include effective CP violation we have obtained results 
for the Chern-Simons asymmetry as a function of CP-violating parameter k. This can now 
be used to estimate the value of k needed to reproduce the observed baryon asymmetry 
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riB/n^- Following the usual arguments, this ratio is related to ub/ s, which is approximately 
conserved in time, where s is the entropy density. After BBN, s = 7.04 I^Oj, whereas 
after the tachyonic electroweak transition it is given by s = (27r^/45) g^T^, = 86.25 the 
effective number of d.o.f. in quarks, leptons, gluons and photons, and T = 4:.Q^Jm}J{GeV) 



GeV from ( 5. 231 ). We also need mw = 80.5 GeV 0. Our results (5/7,5^) for the asymmetry 
{k = Wtt^kiti^^) now lead to 



— = (0.79 ±0.31) X lO^^Km^ {mH = mw), 



(0.46 ± 0.21) X IQ-^K {niH = V2 



mw ■ 



(6.1) 
(6.2) 



Given our approximations of an instantaneous quench, initial temperature zero and the 
neglect of damping effects from degrees of freedom left out in the simulation, we con- 
sider these upper limits on the generated baryon asymmetry^ (the value ( |6.2| ) is proba- 
bly somewhat low, cf. sect. 5.3). To reproduce the observed baryon asymmetry ( [1.1| ) for 
niH = V^mw — 114 GeV, k has to be 



1.4 X 10" 



K ~ 



m 



w 



2.2 X 10" 
ITeV^ 



(6.3) 



which does not seem to be particularly large. Phrasing it differently, venturing the scale 



M = mw (cf. (1-6)), it means 



5cp ^ 0.7 X 10" 



(6.4) 



which is smaller than the Standard Model J ~ 3 x 10 ^ (cf. (1.4)). We see this as encour- 
agement for further pursuing scenarios for ESM electroweak baryogenesis. 
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A. Lattice formulation 

This appendix contains details of the lattice formulation. The lattice spacings in the three 
spatial directions are ai = 02 = as = o, and in the time direction it is oq = at (all are 
positive). The parallel transporter, also called the link variable, from lattice point x + a^/i 
to X is^ C/^,x = u'^_^x+a^fi ~ 6xp[— ia^A|^(a;)r''/2], such that e.g. the forward covariant 
derivative on the Higgs doublet is 

D^(j){x) = [U^^x(t>{x + a^fi) - (t){x)]/a^. (A.l) 

^Note, however, that we have also not taken into account the effect of CP violation before the transition. 
This may bias the initial distribution into an asymmetric form that still has zero average baryon number. 
See also [||, g. 

^In this appendix we use summation convention for the group indices a, fo, . . ., but not for the space-time 
indices ^, m, n, . . .. 
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In the following we use a re-scaled matrix form of the Higgs field defined by 



-(t>*i{x) ^2{X) 



, ^Tr = Xa^ ^t^ 



Furthermore, we shall use lattice units for x and the covariant derivatives, e., 



1. 



(A.2) 



(A.3) 



The backward covariant derivative (indicated with a prime) reads 
A.l Action 

We discretize the action on a space-time lattice of size L^t = (Na)'^ x Ntat with at/a ^ 1 



(A.4) 



in the following way |35|: 



X \_ n ^ ' n>m ^ ' 

-f P'h^Tt: (Do$x)^ (i^o^aO -pfiY^lTr (Dn^x)^ {Dn^x) 



( ^Tr [cDt - ^li,") - /3,iTV $,]Tr [FFJ^,. 



where C/on,x and Umn,x are timelike and spacelike plaquette fields defined as 

Ufj,u,x = U},^.j. = U ^^xUu,x+jjU'lj_^x+u^i',x- 

The CP-violating term contains 



where Unu,x is a symmetrized plaquette field: 



Matching continuum and lattice parameters we find: 

Rt _ ^ <^ as _ ^ at _ 1 a as _^at 
g at g a A at a a 



A a 



R — ^ 2 



{amuf 



(A.5) 



(A.6) 



(A.7) 



(A.8) 

(A.9) 
(A.IO) 
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A. 2 Equations of motion 

Varying the action with respect to A'^ ,^. gives the classical equations of motion: 



Ph 



2^2 



am 



H 



Ph 

Jj/abrj. \. brj 



/3f 



G 



t 

n,a:+0 



Tr 



a;+n+0 a:+n+0 



(A.11) 



+ ( Tr 



xTr 



)} 



(A.12) 



Above, we have defined the lattice E-field ^ through 



.X = iTr 



(A.13) 



We have chosen the temporal gauge, Uq^x = 0, and D'^ is the backward covariant derivative 
in the adjoint representation, e.g. 



7->/a6 Tpb rpa 



1. 



-Tr 



'~^-m,x' '^-m,xi 



■m ^j,x ^j,x 2 

The variation with respect to Aq gives the Gauss constraint: 



j,x—m' 



(A.14) 



G 



(A.15) 



2^ 



A Noether argument tells us that the quantity is conserved in (lattice) time for each x 
when evolving the system with the equations of motion. 

Because of the symmetrization in the discretization of Tr FF, which we found to be 
important for reducing discretization errors, both the Higgs and the gauge field equations of 
motion are implicit. They depend in a non-linear way on forward (in time) links and Higgs 
fields. We solve these equations by iteration. With sufficiently small time steps and/or k 
convergence should be good, and indeed we encountered no problems, iterating between 
3 and 9 times per time step with at /a = 0.05. However, the computational costs were 
correspondingly large. We iterated to computer accuracy (double precision) and checked 
that the Gauss constraint is satisfied, again to computer accuracy. 
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A. 3 Chern-Simons number and winding number 

The lattice expression for tlie change in Chern-Simons number is given by 

1 * 

Ncsit) - NcsiO) = r^Y^Y.^ [FFU,., (A.16) 

^ X 

where Tr [FF]iat,x is given in ( |A.7] ). The following lattice implementation of the winding 
number in the Higgs field, 

where Vx = V{x) with V{x) the SU(2) matrix defined in (|5.5| ), turned out to perform 
satisfactorily in our simulation. 

B. Coping with the Gauss constraints 
B.l Initial Higgs fields 

Our task is to generate an ensemble of initial conditions for the Higgs field in our classical 



simulations according to the distribution (^.4.5; 4.14| ). The distribution applies to a real 
scalar field with the definitions 



k 



^0k, vr(x) = j:^vr,. (B.l) 



In finite volume the momentum label takes on the discrete values k = 27rn/L, and on the 
lattice with lattice spacing a = L/N we can choose Uj = —{N/2 — 1),...,0,..., N/2. The 
reality of the fields imposes (f)]^ = 4>^_^ and similarly for the canonical momenta, except for 
the special modes n = (0 or N/2, or N/2, or N/2), for which exp(ik • x) is real, and 
so the reality conditions imply that the corresponding 0k and vTk are real. We shall term 
them "corner" modes. There are 8 of them. The rest we may as well name "bulk" modes. 
We can write these latter complex variables as 

^ K + K), < = a/^(c£ + 0- (b.2) 



V^^k V ^ 

We have put in a label a = 1,2,3,4 for the real fields of the Higgs field doublet: 



V2 ' V2 

and similar for the canonical conjugate vr. On the lattice we have 



(B.3) 



ujk ^ + kf^,, kL = ^ (2 - 2cos(/c,a)) /a\ (B.4) 

3 

and the lattice Higgs doublet is defined as 



(t)^ = a(t){x). (B.5) 
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We continue to use lattice units and define the lattice canonical conjugate of (/>x as 



5o0x = (A 



x+0 



(B.6) 



For the corner modes, the and are just zero. 

For each spatial lattice site there are three Gauss constraints (one for each generator 
of the gauge group) to be satisfied by the initial configurations. In terms of'^ 



TV [ao$V°$]x, a = l, 



G 



the constraints read 



n.x 



(B.7) 



(B.8) 



where we have used the fact that our initial A'^ = 0. Integrating the left and right hand 
sides in a periodic volume, this reduces to three global Gauss constraints, 



,/3 



(3 = 1,2,3. 



(B.9) 



These are three constraints on our many Fourier modes. Because we have four real scalar 
fields it is convenient to impose and additional fourth constraint: 







E/'x, /? = 4, 



(B.IO) 



which can be seen as the requirement that the total hypercharge also vanishes in the 
periodic volume. Since it is only one extra global constraint on the many modes we do not 
expect that imposing it or not has a noticeable effect on the physical results. 

Inserting equations ( |B.1| ), (|B.2| ) and (B^) into ( |B.9| ) and ( B.IO ) we get a set of con- 
straints on the Ok, 6k) Ck, dkj of which the last one (/3 = 4) is 



kbulk 



2(« + «) 



where we have singled out the k = (0, 0, 0) component for later use, and dumped the rest 
into the quantity f^. In a similar way it is easy to find for the constraints for /3 = 1,2, 3: 



= cofflo - clal + cofflo - CqOo + fi, 
= Cofflo - Coao + Cofflo - CqOo + f2, 
= CQflo - 4^0 + 4^0 - ctah + h, 



(B.12) 
(B.13) 
(B.14) 



where again /i,2,3 are straightforward but long combinations of the remaining variables. 
These equations are linear in all the variables and can be written as 



= Mco + f , 



(B.15) 



^We denote the 2x2 unit matrix by r* 
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with M = M(ao) a 4 by 4 matrix. We are left with the task of generating sets of random 
numbers according to the distribution 

_ i ^ "k + \ + + 5(Mco + f). (B.16) 



We can integrate out the Cq to get 



P(a,b,c,d) (X - — rrexp — > — r 1- \ , (B.17) 

^ ' ' ' ^ detM ^ I 2 ^ nfc + 1/2 2 no + 1/2 1' ^ ^ 

where the prime indicates that the sum over modes no longer includes the Cq, and 

= {M-^ffM-^f = (B.18) 
detM = ^ag2^ (B.19) 

a 

The distribution is no longer a product of Gaussians and we sample it using a Monte-Carlo 
algorithm. For every realization of the remaining unconstrained variables, we then solve 
for the Cq, such that the global Gauss constraints are fulfilled. 

Sampling a distribution with a 5-function can be difficult, and can require very elab- 
orate Monte-Carlo algorithms. However, having sufficiently many degrees of freedom to 
"distribute the constraints on" , the deviation from (in this case) a Gaussian distribution is 
small. In our simulation, where there are about 50 unstable modes (when using Just-the- 
half) and the whole collection of lattice modes when using Thermal, we did not encounter 
such problems. 

B.2 Initial gauge fields 

The Gauss constraint (eqs. A.15| ) imposes a constraint on the gauge fields at each point in 



space, given a Higgs field configuration. Keeping in mind that the gauge fields only grow 
large and classical because of the coupling to the Higgs fields, we set the gauge fields to 
zero initially. This simplifies the constraints dramatically, to read: 

= C," = -d'^E^^^, + ^Tr [(9o^x)^ ir^^.] ■ (B.20) 

We put to zero at each point in space, which means that there are no sources other 
than the Higgs field. Given the Higgs fields we can now write: 

dnK,Ax) = P'^, a = 1,2, 3, (B.21) 

which is just a Coulomb equation, and is solved by introducing a set of scalar functions Xx 
so that 

= -dnXl- (B.22) 
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Fourier transforming 

-dnd'^xl = Px (B.23) 

gives 

r(k) = ^. (B.24) 

^lat 

Since we use a periodic volume, should vanish at k = 0, i.e. the global Gauss constraints 
on the Higgs fields 

Q = Y^dX^^ = Y,Pl (B.25) 

X X 

have to be satisfied. Their implementation is described in the previous section. 
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